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An interesting characteristic of magnetospheric chorus is the presence of 
a frequency gap at cj ~ O.Sffg, where ffe is the electron cyclotron angu¬ 
lar frequency. Recent chorus observations sometimes show additional gaps 
near O.Sfle and 0.6fie- Here we present a novel nonlinear mechanism for the 
formation of these gaps using Hamiltonian theory and test-particle simula¬ 
tions in a homogeneous, magnetized, collisionless plasma. We hnd that an 
oblique whistler wave with frequency at a fraction of the electron cyclotron 
frequency can resonate with electrons, leading to effective energy exchange 
between the wave and particles. 
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1. Introduction 

Chorus in the Earth’s magnetosphere is a type of whistler wave that has been analysed 
for decades using ground-based and satellite observations. Intensive studies have gone 
into understanding the generation and propagation of these waves, and how they affect 
the magnetospheric plasma [Sazhin and Hayakawa, 1992, and references therein]. Chorus 
is believed to be excited by cyclotron resonance with anisotropic {T± > Ty, where _L 
and II correspond to directions perpendicular and parallel to the background magnetic 
held) electrons with energy > 1 keV injected into the inner magnetosphere. Studies have 
shown that the resonance between chorus and relativistic electrons plays an important 
role in radiation-belt dynamics [Thorne, 2010, and references therein]. For example, pitch 
angle scattering by chorus is a major loss mechanism for trapped electrons in the outer 
radiation belt. Local acceleration due to interactions between chorus and electrons inside 
the radiation belt may be a major mechanism for enhanced relativistic electron population 
[e.g. Thorne et ah, 2013]. 

An important feature of chorus is the presence of a gap at one-half the electron cyclotron 
frequency in its spectrum, separating two frequency bands (therefore called “banded cho¬ 
rus”), a lower band with 0.1 < u/flg < 0.5 and an upper band with 0.5 < u/Qe < 
0.8 [Meredith et ah, 2012; Li et ah, 2013], where hie is fhe electron cyclotron angular 
frequency. Such banded chorus is not unique to the Earth’s magnetosphere, but is also 
observed in Saturn’s magnetosphere [Hospodarsky et ah, 2008]. Recently, using Clus¬ 
ter spacecraft measurements, Macusova et ah [2014] reported that sometimes chorus can 
have more than two bands, with additional gaps near O.Ohle and O.Ohle- Most of these 
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“multi-banded” chorus events were observed with oblique wave normal angles during dis¬ 
turbed geomagnetic conditions. An example of multi-banded chorus observed by Van 
Allen Probes A on Feb 10, 2013 is shown in Figure 1. 

The narrow gap at O.Shle between the two bands of chorus has received a lot of attention 
since it was discovered (for a review of existing theories, see Sazhin and Hayakawa [1992] 
and Liu et ah [2011]). Liu et al. [2011] hypothesized that two bands of chorus are excited by 
two different electron populations with temperature anisotropies through linear instability, 
and the hypothesis has been tested in a case study using Van Allen Probes data by Fu 
et al. [2014]. Schriver et al. [2010] explored the possibility of generating the lower band 
through nonlinear wave-wave coupling of the upper band chorus. Bell et al. [2009] assumed 
that different bands are generated in ducts of either enhanced or depleted cold plasma 
density. Omura et al. [2009] explained the gap at O.Sfle, where the group velocity of the 
whistler wave equals the phase velocity, as nonlinear damping of a slightly oblique whistler 
wave packet propagating away from the magnetic equator, taking into account the spatial 
inhomogeneity of the magnetic held. Despite intensive research, a conclusive explanation 
for banded chorus has not been given yet. 

Existing models for banded chorus have only considered the primary (linear) resonances 
when the resonant particles described by unperturbed orbits see a “time-independent” 
wave held. However it was recently found that the nonlinear resonances which develop 
by taking into account the perturbed particle motion in the wave held can also prompt 
energy exchange between the wave and particles. For instance, island overlapping due to 
nonlinear resonances can cause stochastic ion heating in an oblique Alfven wave with a 
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sub-cyclotron frequency in the solar corona [Chen et ah, 2001; Guo et ah, 2008]. Can the 
nonlinear resonance happen between chorus and electrons? If it can, how will it affect the 
electrons and chorus? 

Inspired by the interesting feature of nonlinear resonances at sub-cyclotron frequencies, 
we propose the following scenario for the banded chorus observation illustrated in Figure 
1. A warm (a few hundred eV), anisotropic (T± > Ty) electron velocity distribution 
drives the whistler anisotropy instability [Gary et ah, 2000, 2011], which gives rise to 
continuous narrowband (O.dhle < OJ < 0.7fie) enhanced magnetic spectra that appear as a 
relatively coherent temporal waveform. The cold electron (1-100 eV) response, as will be 
shown by both theoretical analysis and test-particle simulations in this work, demonstrates 
that there is a nonlinear wave-particle interaction whereby certain electrons come into 
subharmonic resonance with certain Fourier components. If this interaction transfers 
energy from the fluctuations to the electrons, the resonant Fourier components will be 
damped and the fluctuation spectra will develop gaps at fle/2 and other subharmonics as 
shown in Figure 1. 

In this paper, we show that, in the absence of primary resonances, an oblique whistler 
wave with a frequency at a fraction of fie is able to resonate with the cold electrons 
nonlinearly, leading to nonlinear damping/growth of the wave with certain electron dis¬ 
tributions. This nonlinear mechanism, which involves only wave-particle interactions and 
works in homogeneous plasmas with a uniform magnetic held, can provide a complemen¬ 
tary element to existing theories on chorus. In addition, it can explain additional gaps in 
chorus spectra around O.Shle and O.bhle, as reported recently by Macusova et ah [2014]. 
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The rest of the paper is organized as follows. In Section 2, we present a theoreti¬ 
cal framework for analysing the dynamics of electrons in an oblique whistler wave with 
uniform background magnetic held. The structures of nonlinear resonances are analysed 
using Poincare maps by solving the equations of motion numerically and conhrmed by our 
theoretical calculations employing the Lie perturbation method. In Section 3, we show 
the results of test-particle simulations for an ensemble of electrons with certain velocity 
distributions. The effects of nonlinear resonances on the electron distribution function 
and the total kinetic energy are investigated. In Section 4, we discuss how this nonlinear 
mechanism is related to frequency gaps in magnetospheric chorus. Finally, conclusions 
are given in Section 5. 

2. Theoretical Analysis 

In this section, we analyse the dynamics of electrons in a single oblique whistler wave 
and an uniform background magnetic held using Hamiltonian theory. 

2.1. Hamiltonian 

For simplicity, we consider a uniform plasma in a uniform background magnetic held, 
Bq = Bqt., and the whistler wave dispersion relation in the cold plasma limit is given as 
[Stix, 1992, Equation (2-45)] 



( 1 ) 


where Ue is the plasma frequency, k is the wave number, and the whistler wave is oblique 
with a small wave normal angle 9 with respect to the background magnetic held so that 
a = tan0 = k^/kz. Assuming 9^1 and using the cold plasma theory [Stix, 1992], the 
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vector potential can be written as 

A = Bqxy + e(-Bo/fc^)(sin \l/x + cos 0 cos \l/y), (2) 

where \I/ = kxX+kzZ—ut is the phase of the wave and e = B^/Bq denotes the pertnrbation 
magnitnde. To eliminate electric helds, we move to the wave frame by a transformation, 
x' = X— {u/kz)tz. Normalizing time to l/flg, magnetic held to Bq, mass to nie and length 
to 1/kz, the normalized Hamiltonian for the electron is given by 
H = ^[pI + {Py + xY +pY\ + e [Px^mip + {py + x) cos 6^ cos V'] + ^^^(1 “ sin^ 0 cos^ , 

(3) 

where Py, and Pz are the canonical momenta and ip = ax + z. In the wave frame, the 
electron energy is conserved as H does not depend on time explicitly and Py is a constant 
of motion since H is independent of y. A set of eqnations of motion for the electron can 
be readily obtained from the Hamilton’s eqnation [e.g. Jose and Saletan, 1998]: 


X = Px + esin T, (4) 

Px = ~{Py + x) — e[apx cos T — a(py + x) cos 0 sin T + cos 9 cos T] 

.2 

——a sin^ 6* sin(2T), (5) 

^ = Pz, ( 6 ) 

g2 

Px = —e[pa; cos T — (py + x) cos6* sin T] — — sin^ 6'sin(2\l'), (7) 


where dots represent time derivatives. 

Withont loss of generality, we choose the following parameters relevant to chorns in the 
Earth’s magnetosphere. The whistler wave satishes u = 0.4He, k = 0.9a;e/c and 9 = 26.6° 
(a = tan 6* = 0.5), where c is the speed of light and oJe/kle = 5. So the parallel phase speed 
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of the wave is uj/kz ~ 0.1c, and the energy is normalized to me{^e/kzY ~ O.OOmeC^ = 
30.66 keV. Note that the choice of the whistler wave frequency is for illustration purpose. 
It will be shown later that, for waves closer to O.Sflg, the difference is only in the energy 
of resonant electrons. The extension to a whistler wave packet with multiple Fourier 
components is also straightforward [Lu and Chen, 2009]. Unless otherwise specihed, the 
above parameters are used in the rest of the paper. 

2.2. Poincare maps 

In our Hamiltonian model, electrons are moving in a 4-dimensional phase space 
z,Pz). With Poincare surfaces of section (or maps), one can visualize the wave- 
particle resonances in phase space. We construct a Poincare map in {pz,z) by recording 
points when the trajectory of an electron in phase space crosses the surface of x = 0 with 
Px > 0. Electrons are initialized with x = 0, Py = 0, hxed H and a range of parallel 
velocities Pz^s. In the top panel of Figure 2, a map for electrons with energy H = 0.3 is 
shown in the presence of an oblique whistler wave with e = 0.02. In the map, a major 
island in (p^, z) plane located at = 0 demonstrate trapping of particles by the wave 
satisfying cj — kzVz = 0 (note pz = meivz — oj/kz) before normalization), which is the 
well-known Landau resonance condition. Another set of two islands develop at p^ = — 
which corresponds to a new resonance condition 

u - kzVz = kle/2, ( 8 ) 

suggesting that electrons can resonate with hnite amplitude oblique whistler waves at 
half cyclotron frequency. A similar mechanism has been proposed for the sub-cyclotron 
resonance between ions and oblique Alfven waves [Chen et ah, 2001; Guo et ah, 2008]. 
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With higher energy, the primary cyclotron resonance between electrons and the whistler 
wave can be seen from the island located at Pz = (not shown). 


2.3. Sub-cyclotron resonances 

We can show that the half cyclotron resonance is the result of nonlinear dynamics by 
taking into account the perturbed electron orbit in the presence of wave helds. The anal¬ 
ysis is greatly simplihed in the so-called guiding-center coordinates. Through a canonical 
transformation [Guo et ah, 2008], the Hamiltonian becomes 


H = Ho+ + e^H2 


(9) 


where Ho = J + pl/2 is the unperturbed guiding-center Hamiltonian and 


Hi 

H2 


cos (j) sin(Q!\/^J sin (j) + z) + sincj) cos 9 cos{ay /^sin (j) + z) 
— sin^ 9 cos^(a\/^ sin (j) + z). 


( 10 ) 

( 11 ) 


Here, J = [{py + x)^ -|- p^]/2 denotes the perpendicular energy and (j) the gyro-angle. 
Interestingly, if the whistler wave is purely parallel to the background magnetic held 
(0 = 0), H 2 vanishes and Hi = \/^sin(z — 0) contains only the primary resonance. So 
the electron dynamics becomes integrable and there will not be any nonlinear resonance. 
Due to the dependence of the perturbed Hamiltonian on 0 and 2 ;, both J and pz vary in 
time. These variations are related to the magnetic drift of the guiding center which can 
then give rise to nonlinear resonances via the term with being the magnetic 

drift. 

To analyse the perturbed Hamiltonian, we introduce the Lie transform method [e.g. 
Jose and Saletan, 1998] which utilizes particular canonical transformations to obtain per- 
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turbation series in terms of Poisson brackets. Similar to the approach in Guo et ah [2008], 
a generating function 



cos(n 0 + z) 

n + Pz 


( 12 ) 


is obtained to remove the hrst order perturbations by setting [Wi, Hq] + Hi = 0 with [, ] 
being the Poisson bracket. This generating function indicates singularities when = —n, 
with n being any integer. In the laboratory frame, they simply correspond to the well 
known cyclotron resonance condition u — kzVz — nQe = 0 when n 7 ^ 0 and the Landau 
resonance condition u — k^Vz = 0 when n = 0 . 

After the Lie transformation, the new Hamiltonian becomes H' = Hq + €^{H 2 + 
[Wi, Hi\/2) + O(e^), where the second order perturbation reads 


H' = 


^ sin 0 J^riJn) cos^/(j/> T 2z^ ^ “t“ siu 0 cos(?7r Tl^cj) 


m,n 


m,n 


(13) 


with I = m + n. The functions Am,n{J,Pz), Bm.,n{.JiPz) are introduced for simplicity. Their 
explicit expressions are tediously long and will be presented later in a separate paper. 
The new Hamiltonian clearly shows the existence of nonlinear resonances when the phase 
{l(j) + 2z) remains a constant or pz = —1/2. In the laboratory frame, they correspond to the 
second order resonance conditions uj — k^Vz — {l/2)Qe = 0. The second term in Equation 
(13) does not contribute to resonances. Therefore, our analytical calculation conhrms the 
existence of the nonlinear resonance at the half cyclotron frequency (/ = 1). Following 
the result in Guo et ah [2008], we know that the island width of second order resonances 
is proportional to e while it is y/e for the primary resonances. This means the damping 
due to second order resonance is weaker than the primary resonance by square root of the 
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wave magnitude. Continuing to the third order expansion, we will obtain resonances at 
Pz = —1/3 which is present in the bottom panel of Figure 2 where the wave has a larger 
wave amplitude. 

3. Test particle simulations 

A wave-particle resonance causes efficient energy exchange between the wave and reso¬ 
nant particles. For example, a wave is damped via Landau resonance while the resonant 
particles gain same amount of kinetic energy the wave loses. To further illustrate how the 
nonlinear resonance, at ca — — VL^/2 = 0, can damp waves around the half cyclotron 

frequency, we perform test-particle simulations of a large number of electrons and inves¬ 
tigate the changes in the electron distribution and their kinetic energy in the presence of 
an oblique whistler wave. 

In our test-particle simulation, about 10® electrons with Maxwellian velocity distribution 
in the lab frame. 



(14) 


are loaded initially. Here we choose electron temperatures = Ty ~ 300 eV so that 
Vt = y/T\Jrn/ |(i^ — f^e/2)/fc|||. Electrons are advanced according to equations of 
motion given by Equations. (4)-(7). In the presence of an oblique whistler wave with wave 
amplitude B^/Bo = 0.02, the energy (calculated in the lab frame) evolution of these 
electrons is shown in Figure 3a. After initial oscillations, the averaged energy per electron 
(E) increases at a rate AE/^EqAI) k, 10“^He (where Eq is the averaged initial energy 
per electron) and the energy gain is predominantly in the perpendicular direction. In a 
self-consistent simulation where the wave held follows Maxwell equations and the total 
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energy of the system is conserved, the increase of the electron energy must come from the 
decrease of wave energy. This implies that the wave will be damped by the nonlinearly 
resonant electrons. The rate of particle energy gain normalized to the wave energy is 
2fioneme{^le/kzyAE/{BlAt) ^ 0.02hle- 

The electron velocity distributions in the wave frame at tQe = 0, 60,120,180 are shown 
in Figure 3b. Clearly, the parallel velocity distribution f{pz) deviates from the initial 
Maxwellian in the vicinity of pz = —0.5, which corresponds to the second order resonance 
shown in the Poincare map in Figure 2. The changes in both f{px) and f{py) are small (not 
shown). For the given Maxwellian distribution function, the nonlinear cyclotron damping 
steepens the gradient of f{pz) around pz = —0.5, causing the parallel kinetic energy to 
decrease, and the perpendicular kinetic energy to increase (see Figure 3b and 3c). This 
is similar to the primary cyclotron resonance where the electrons are scattered along the 
constant energy surfaces in the wave frame as pointed out by Kennel and Engelmann 
[1966]. 

Just like the Landau or cyclotron resonance, the nonlinear resonance can lead to damp¬ 
ing/ growth of the waves at sub-cyclotron frequencies when the electron distribution satis- 
hes certain conditions (the stability criteria). Even though the particle motion is described 
as periodic in our theoretical analysis and test-particle simulations, what determines the 
wave damping/growth rate is how fast the resonant particles can take (give) energy from 
(to) the perturbed helds. This is the same as in Landau or cyclotron damping, which 
requires the resonant particles to gain energy at a rate faster than their bounce frequency 
in the perturbed held. The wave begins to saturate when this condition is violated. The 
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difference in the nonlinear damping mechanism is that the relevant bonnce motion is now 
described by the second order pertnrbed Hamiltonian and thns becomes mnch slower than 
the linear bonnce motion. A test-particle simnlation describes the particle dynamics at 
the beginning stage of the wave-particle interaction and can predict the rate of energy gain 
or loss by resonant particles depending on the distribntions, as discussed in the following. 

When the temperature anisotropy T±/T\\ is strong enough, the nonlinear resonance may 
also lead to decrease of the particle energy and thus instability of the whistler wave. This 
is demonstrated by a comparison of runs with different initial anisotropies, as shown in 
Figure 4. Same parameters as those in Figure 3 except the perpendicular temperature are 
used in these runs. For Tj_/T|| = 100, the electrons lose energy, which implies wave growth 
in a self-consistent situation. Interestingly, for anisotropies well below the instability 
threshold, the temperature anisotropy can enhance the nonlinear wave-particle resonance, 
as shown by the case with = 3 initially. This is consistent with our analysis as in 

Equation (12) where a factor of \/^ shows the dependence on the particle perpendicular 
energy. 

To further illustrate the sub-cyclotron resonance and clarify that the energy gain of 
electrons is not dominated by numerical heating or the pseudo-heating [Dong and Singh, 
2013] in our simulations, we carry out two comparison test cases. In the hrst case, we 
artihcially load electrons with f(pz) hat in the range —0.6 < pz < —0.4 as shown in 
Figure 5(a). Due to the nonlinear resonance near half cyclotron frequency, electrons 
gain energy from the wave. In the second case, we remove the resonant particles within 
—0.55 < Pz < —0.45 from the distribution as shown in Figure 5(b). The resulting energy 
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gain in the lab frame (Fignre 5c) is significantly rednced. Therefore it is clear that the 
dominant energy gain of electrons is caused by the nonlinear resonance. 

4. Application to Banded Chorus 

In the Earth’s magnetosphere, cold plasmaspheric electrons (~ 1 eV) extend into a 
region called the “plasma trough” with L = 4—8 (L is the equatorial distance of a magnetic 
held line from the center of the Earth in the unit of Earth’s radius) [e.g. Carpenter and 
Anderson, 1992], where they meet ~ 1 keV electrons from the plasma sheet. Banded 
chorus is frequently observed in this region [e.g. Meredith et ah, 2012]. We conjecture 
that the gap of the banded chorus may be caused by the nonlinear damping of oblique 
whistler waves by these cold electrons with energy 1-100 eV. For 10 eV electrons, whose 
thermal velocity Uf ~ 6 x 10“^c, the nonlinear resonance condition (Equation 8) becomes 
uj ^ O.Shle, i-e., these electrons could damp oblique whistler waves with frequency near 
O.Srie, leaving a gap in the wave power spectrum. 

Wave analysis often shows chorus is quasi-parallel, but we have also shown that even 
with a wave normal angle of 26.6°, nonlinear wave-particle interactions can be signihcant. 
Furthermore, even if the chorus is excited purely parallel in the source region close to the 
magnetic equator, as the wave propagates away from the equator, it will soon have an 
oblique component due to the curved nature of Earth’s dipole-like magnetic held. 

Here, we have presented cases with a modest wave amplitude B^j/Bq = 0.02 to better 
illustrate the nonlinear resonance mechanism with less computational constrains. The 
numerical integration scheme can introduce errors into the particle’s energy, therefore, for 
smaller wave amplitudes, we need to both reduce the time step to suppress the numerical 
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heating and increase number of particles to resolve the resonance structure in phase space. 
Typical amplitude of chorus in the magnetosphere is about B^j/Bq ~ 0.001, but large 
amplitude whistler waves have also been reported with B^/Bq > 0.01 [e.g. Santolik et ah, 
2014]. The sub-cyclotron resonance can be strong and effective in damping these large 
amplitude whistler waves as shown in the previous section. For smaller wave amplitude, 
we have also carried out test-particle simulations with B^/Bq = 0.002, uj = O.dOfle; ^|| = 
50eV, and T± = lOOeV. The results (not shown here) indicate that the same nonlinear 
resonance occurs, although the rate of the electron energy gain is reduced {AE / (EAt) ^ 
6 X 10“®r2e) due to the lower wave amplitude. Therefore the proposed nonlinear damping 
mechanism at sub-cyclotron frequencies is robust and can play a role for waves observed in 
the magnetosphere. Furthermore, in the presence of multiple whistler modes (or a narrow 
band) near half cyclotron frequency, the nonlinear damping is expected to be enhanced 
due to the overlapping of resonances as suggested by Chen et ah [2001] and Lu and Chen 
[2009]. This effect will be addressed in a separate paper. 

If the wave amplitude increases to B^/Bq = 0.05, additional resonant islands near = 
— 1/3 and pz = —2/3 develop in the Poincare map (Figure 2, bottom panel), which may 
explain the additional gaps near O.dfle uud O.Ohlg observed by Macusova et ah [2014]. As 
Macusova et ah [2014] reported, multi-banded chorus were observed during more disturbed 
times with the average Kp ~ 3, larger than the average Kp ~ 2 for chorus observed with 
one or two bands. These higher order nonlinear resonances, as well as the half cyclotron 
resonance we have shown, belong to a set of nonlinear resonant conditions below the 
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cyclotron frequency, cu — k\\v\\ = NQ^/M, where M, N are integer and N < M [Guo et al., 

2008], 

5. Conclusions 

In this paper, we have presented a nonlinear resonant mechanism between an oblique 
whistler wave and electrons, satisfying the resonant condition ca — k\\v\\ = Ge/2, by the¬ 
oretical analysis and test-particle simulations. Our mechanism works in a homogeneous 
plasma with a constant background magnetic held, and may explain the frequency gap at 
O.SGe frequently observed in the power spectra of magnetospheric chorus. Further more, 
similar nonlinear resonances may explain the frequency gaps at O.SGe and 0.6Ge in chorus 
observations recently reported by Macusova et al. [2014] and as shown in Figure 1. This 
mechanism provides a complementary element to existing theories on banded chorus. The 
detailed theoretical analysis (Section 2) will be presented in a separate paper later. The 
ability of this nonlinear mechanism to explain frequency gaps in chorus emissions needs 
to be further investigated in a self-consistent way (e.g. via particle-in-cell simulations) to 
address the dependence of damping rates on various plasma and wave parameters. 
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seconds 


Figure 1. A multi-banded chorus measured by Van Allen Probes A on Feb 10, 2013: the 
magnetic held spectrogram with two frequency gaps at 0.5/ce and 0.6/ce, where fee is the local 
electron gyrofrequency. Black dashed lines indicate 0.4/ce, 0.5/ce and 0.6/ce, respectively. 
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Figure 2. (top) Poincare map in {pz,z) plane for electrons in the presence of an oblique 
whistler wave with an amplitude B^/Bq = 0.02. The energy of electrons in the wave frame 
is hxed with H = 0.3. Islands develop near Pz = 0 and p^ = —(bottom) Similar Poincare 
map, but with a larger wave amplitude B^/Bq = 0.05. Additional resonant islands develop near 
p^ = and p^ = 
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B„IBo = 0.02 




Pz 

Figure 3. (a) Change of the averaged energies per test electron in the presence of an obliqne 

whistler wave with an amplitnde B^/Bq = 0.02. The averaged initial total energy per electron 
.E'tot,o = 0.028. (b) The distributions of f{pz) at different times of the simulation. Changes of 
f{px) and f{py) are minor. 
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Figure 4. Relative change of energies per electron for cases with temperature anisotropy 
T}_/T\\ = 1, 3,10,100. All parameters except are fixed. 
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Figure 5. The evolution of the distribution function f{pz) for (a) an initial flat-top distribution, 
and (b) a flat-top distribution with resonant particles removed. The comparison of the change 
of total energy per electron is shown in (c). 
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